Exact cosmological solutions of models with an interacting dark sector 



o 

(N 
+-> 

o 

O 

m 



A. B. Pavan °E Elisa G. M. Ferreira 6 S. Micheletti 6 1 J. C. C. dc Souza c l and E. Abdalla^ 

a Institute/ de Ciencias Exatas, Universidade Federal de Itajubd, 
Av. BPS 1303 Pinheirinho, 37,500-903, Itajubd, MG, Brazil 
b Instituto de Fisica, Universidade de Sao Paulo, 
CP. 66318, 05315-970, Sao Paulo, SP, Brazil and 
c Centro de Ciencias Naturais e Humanas, Universidade Federal do ABC, 
Rua Santa Adelia 166, 09210-170, Santo Andre, SP, Brazil 
(Dated: October 25, 2012) 

In this work we extend the First Order Formalism for cosmological models that present an interac- 
tion between a fermionic and a scalar field. Cosmological exact solutions describing universes filled 
with interacting dark energy and dark matter have been obtained. Viable cosmological solutions 
with an early period of decelerated expansion followed by late acceleration have been found, notably 
one which presents a dark matter component dominating in the past and a dark energy component 
dominating in the future. In another one, the dark energy alone is the responsible for both periods, 
similar to a Chaplygin gas case. Exclusively accelerating solutions have also been obtained. 
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I. INTRODUCTION 



The recent improvements in the observational techniques available for the measurement of cosmological parameters 
show, with increasing accuracy, that the universe is composed mainly of two mysterious entities, the so-called dark 
(3J^ energy and dark matter. The first one is believed to be the cause of the observed accelerated expansion of the universe 
[2-01 1 and corresponds to approximately 70% of its total energy density. The latter corresponds to almost 25% of 
the energy density of the universe and plays an important role in large structure formation. The reader is referred to 
[H-[IH for a review on the theoretical developments that followed these observations. 
" An interaction between these two dark sectors is rather plausible, and can be even considered to be a necessary 
^vq , feature in cosmological models based on Quantum Field Theory. An extensive literature in the last years has shown 
■ both theoretical and phenomenological aspects of the coupling between dark matter and dark energy (see, e.g., fltl ]- 
^sO ! [H| ) . The present work is an attempt to solve the equations of motion for a theoretical model based on an interaction 
Lagrangian, in which dark energy is represented by a canonical scalar field, and dark matter takes the form of a 
fermionic field. 

Solving these equations, even for the case in which no interaction is taken into account, is very often an extremely 
arduous task, if not totally impossible, without the use of certain approximations. Many papers rely on intensive 
computer simulations to gather some insight about cosmological parameters, while others tackle the difficulties by 
analysing the phase space dynamics represented by the generic equations of motion for the cosmological model in 
! question. 

In all cases, the viable cosmological solution including dark energy must obey that the energy density of scalar field 
must remain subdominant during radiation and matter dominant eras to allow for structure formation, becoming 
dominant in late times accounting for the current acceleration of the universe. This can be obtained by scaling 
solutions where the energy density of dark energy mimics the background fluid (radiation or matter) [3?J ■ The field 
must also exit this scaling solution in order to describe the late accelerated solution. Models with interaction between 
the dark sectors can account for the exit of scaling, by dynamically modifying the scalar field potential, and satisfy 
all the requirements to describe a viable dark energy scenario [38l . |39| . 

Another method used to obtain exact solutions in cosmology is the First Order Formalism (FOF), developed by 
Bazeia et. al. [4(| , where cosmological solutions with dark energy modeled by a single scalar field were obtained exactly. 
The central point of that formalism is to assume the Hubblc's factor to be a function of the scalar field (j>, thus reducing 
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the second order equations to first order ones. In another work Bazeia et. al. [4l| generalized the first order formalism 
describing a universe filled with dust and dark energy. It is interesting to point that, in their work, a coupling between 
dark matter and dark energy arised naturally. The authors argued that this interaction is an effect of the applied first 
order formalism. 

The aim of this work is to apply the FOF for a model where a interaction between the dark sectors is present. We 
consider a universe filled only with dark matter and dark energy which means that our models only aims to model 
the late universe evolution. We have generalized Bazeias's cosmological models by explicitly adding a fermionic field 
ip, which plays the role of dark matter, and an interaction term between the fermionic and the scalar parts. The 
imposition of such interaction renders a more involved Fricdmann equations and equations of motion, increasing the 
number of variables of the problem and reducing the constraints of the system. Since we have this extra freedom in 
the model, we impose some restrictions to the scale factor, via an Ansatz that relates it to the scalar field itself. 

Even though these restrictions may reduce the possibilities of expansion of the universe, we show in this paper 
that we can still construct a large class of exact solutions by applying the FOF for the interacting models with the 
restriction of the solution for the scalar field as a function of time, <fi(t), to be invertible. This aspect is actually 
essential for solving the equations exactly and is taken for granted in much of the literature on this subject. A few 
works involving dynamical analysis discard this requirement, but they are able to find exact solutions after making 
some different restrictions and approximations based on the structure of the phase space for the models [37| . 

The solutions obtained present only accelerated, and decelerated and accelerated periods. This is in accordance 
with what happens with the phase space analysis, where in most cases we obtain non-viable cosmological solutions 
and the scaling solution, that describes the right cosmological evolution, is harder to find, obtaining only a few for 
each model. We also note that, together with the dark sectors interaction proposed as a hypothesis, the method 



introduces an interaction, as pointed in [41( . This happens because of the choice that the Hubble parameter and scale 
factor depend only on the scalar field, what renders an effective interaction between dark energy and dark matter. 
We also account for this interaction and verify how this influences the solutions. 

The plan of this paper is as follows: in Section [XT] we introduce our model of interacting dark energy and dark 
matter introducing a fairly general interaction form. In Section IlIII we present the first order formalism and present 
the Ansatz necessary to reduce the order of the differential equations. In section IIVI we show the exact solutions 
obtained for different parameters of the formalism. The last section shows our final remarks. 

II. INTERACTING DARK ENERGY AND DARK MATTER MODEL 

First let us show the model we are going to study. We use the Friedmann-Lemaitre-Robertson- Walker (FLRW) 
metric with null curvature 

ds 2 = dt 2 - a(t) 2 (dr 2 + r 2 dd 2 + r 2 sin(6) 2 d(f> 2 ), (1) 
where a(t) is the scale factor. The action that describes this geometry and the material content is 

,1„ ,„., , „, i 



S = j d 4 x^ + l-d^4> -V((fy+ l - V„V - ffirfl] - Mjnl> + PF(<f>)Wj , (2) 

where <j> is the scalar field and V{<j>) its potential, ip is the fermionic field and F{4>) is an interaction. The constants M 
and /3 are the mass of the fermionic field and the coupling constant, respectively and t is a constant that can assume 
the values ±1, if one wants the action to also accomodate the possibility of a phantom field. The functions F{<j>) and 
V{4>) are not fixed in the beginning. 

We have to solve Friedmann's equations, the equation of motion for the scalar field and Dirac equation. We know 
that the homogeneous Dirac equation in a homogeneous and isotropic spacetime has a simple solution in terms of the 
scale factor and we can see that ipip = a/a 3 , where a is a constant related to the energy density of the dark matter at 
the present phase of the history of the universe (j^. Therefore, the system of equations that remains to be solved is 

<P + 3H(j)+ — = — - — , (3) 
H ' 2 = i^{^+ V ^) + [M-pF{<m^}, (4) 
Mi 




-^{4 2 + [M-PF(<f>Mi>}, (5) 
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where H — - is Hubble's factor and M p the Planck's mass. A prime represents derivative with respect to the scalar 
field <t>. 

We also evaluate the energy-momentum tensor. We evaluate it separately for each component, the scalar field and 
the fermionic field, and for the interaction term [52j . We then obtain the energy density, pressure and equation of 
state in the FLRW space-time. As the interaction term, Ci n t = /3F((f))'ipip, depends on both the scalar field and the 
fermionic field we choose not to include the interaction as a part of any of those fluids. Instead, we separate it as a 
different "component" so we can see how it evolve and how dominant this term is, since F((f>) is not known and the 
method can introduce new interactions. 

The energy density /?, the pressure p and the equation of state oj for the scalar field are given by: 

Pde = -z- + V{<t>) , p DE = — V((p) , w DE = - , (6) 

Zip 

where for the standard field case I = 1 and for the phantom field case I = — 1, leading to a change in the sign of the 
kinetic term in the energy density and pressure and resulting in an equation of state parameter such that w ^ — 1. 
For dark matter we have 



a a —fi4> 

PDM = M — , PDA/ = -P<P^T , WflM = — — . (7) 

a 6 a 6 M 



We can see that in the case with interaction the pressure of dark matter does not vanish for all times but it depends 
on both fields and the interaction constant. This happens because when we have interaction the components do not 
conserve alone, but the sum of all components is what is conserved. However, we can stress that depending on the 
solutions, pdm and ojdm can vanish for a certain time. It is important to notice that this choice of separation of the 
components breaks the standard distinction between the dark matter and dark energy, since now both components 
can evolve as matter, dark energy or both, depending on the solution obtained. This is fundamental to understand 
the mixed role played by each component in the solutions obtained in Section IIVI 

For the interaction " component" , the energy density, pressure and equation of state are given by 

(X Oi 

Pint = , Pint = j#0-o j Wint = -1 ■ (8) 

a 6 a 6 

The constant equation of state equal to —1 shows us that the interaction "component" also makes the role of acceler- 
ating the universe, as the dark energy does. Therefore, if the interaction is too strong in comparison with the other 
components, this would act to accelerate the universe. 

We can also define, as usually, the density parameter 51 = p/p cr it, where p cr u = 3HqT/2 with V = M£/4tt. 

The above system of equations ([3]l5|) describes the dynamics of the variables of our interacting system. We must 
solve this system in order to obtain the evolution of our components. We can see that the interaction term brings an 
explicit dependence on the scale factor in the equations adding a new dynamical variable in the differential equations. 
We can notice, for example, that the Fricdmann equation ((4]) can no longer be used as a constraint equation, as it is 
in the case without interaction. 

Because of these differences the system is more complicated to be analysed. If one wants to compute the solutions 
numerically one needs to perform a phase space analysis, that can only be done for a power law-like potential in 
the case of Yukawa couplings (48j , requiring the imposition of assumptions in the variables to be studied for general 
potentials. Because of these difficulties, we explore a different route to study this model based on [13,|4l|, where we 
rely on important simplifications to obtain exact solutions. We want to check if this method and these simplifications 
can be used and give good results in the context of interacting models. 



III. FIRST ORDER FORMALISM - FOF 



Now, we are going to introduce the formalism that we are going to use. Following the central idea of FOF we make 
the simplifying assumption that: 

H(t)=wm)) ==► H = W*4>, (9) 

where the dot stands for time derivative and W$ = dW/d<f>. We require that 4> must be invertible, so solutions will 
only be possible through this method if the scalar field and its time derivatives are smooth, monotonic functions {42j . 
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Other approaches can present solutions that arc not invertible (see, e.g., [Illllj]), but they are usually valid only 
for very specific regions of the phase space of the cosmological models. This invertibility requirement restricts the 
solutions that can be obtained by the method. 

With this in mind we rewrite the Eqs.([5][5|) substituting H by W{4>). The Eq. ([5]) becomes 

-W+ ^r = # 2 + [M-/3F^)]^, (10) 

with r = ^2- . By inspection of the Eq. © , we see that it is necessary to impose some restrictions on the form of the 
scale factor. In order to solve this equation in terms of <j), we make an extra assumption and choose to rewrite the 
scale factor using the Ansatz 

a(t)- 3 = a<j> n J(4>), (11) 

where a is a real constant, n is an integer and J((j>) an arbitrary function of the scalar field. This expression has a 
general form in a way that allow us to obtain a large class of exact solutions with interacting dark energy and dark 
matter, by choosing convenient integers n and functions J((f>) that reduce the order of the equations of motion. A 
more direct approach, without assuming this ansatz, can also be done as shown in Appendix 1. 
Thus, substituting (fTTj) in (fpQ|) we obtain 

M+w<n _ (12) 



[M - pF((p)}aa J(<f>) 

The potential V n {(j)) associated to the scalar field is calculated using (j4}, resulting in 

3TW 2 6 2 

Vntt) = ~ £ Y~ [M ~ (13) 

We can solve Eq. (fT2"|) as an algebraic equation for <p for each value of n. Some of the roots of this equation can 
be imaginary, and will be discarded. We also restrict the values for n in our calculation to n = 0, 1, 2, .... Such a 
procedure reduces the order of the equations that need to be solved and transforms the equation of motion for the 
scalar field into a constraint equation. The functions W(4>), J{4>) and F(4>) need to satisfy the constraint equation: 

3W j, + 3TWW * _ EM ~ PEMMi 



4> n J{4>)\ = o , (14) 

where we have used Eqs (|12ll3[) and the relation <fi = <f> d<fi/d(j) valid whenever the function <p(t) is invertible. This is 
the point where it is crucial for this function to be a one to one map, what is not always true. 

As we can see, the power of this method is that you obtain two important quantities. The scale factor that is exactly 
what it is not known and we would like to obtain (to verify, for example, that a scalar field with a potential yields 
acceleration) and not to choose (as it is done in the direct approach) ; and the form of the interaction, parametrized 
by F((f>)[53| , the most important feature of the interacting models that is in principle unknown. You only impose the 
form of the Hubble parameter and that the scale factor depends on the scalar field given the form of this dependency. 
The method then provides you the form of the evolution of the system as well as the interaction between dark matter 
and dar energy 

So, the setup to obtain the solution is the following. Given W{<f) and J(4>) and chosen a value for n, we write the 
potential ([!]) with the chosen form for these functions and solve the constraint equation ([5]) to obtain the form of the 
coupling F(<f>). We plug this in (|3|) and solve this algebraic equation with the order given by the value chosen for n. 
With this it is possible to obtain <f>(t) and a(t) and evaluate the energy density and pressure of both components. 

As an example, we can show the form of the equations for n = 0. The first order equation will be given by 

W T JiW+T) 2 - U[M - pF{<l>)] a oJ{<l>) 

(p = -^r ± 2i (15) 



with the potential Vq being 



3TW 2 6 2 

VoW = [M - 0F(<j>)]aaJ(cl>) . (16) 
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The constraint equation relating the functions J, W and the interaction F is given by 

with the scale factor being given by a(t)~ 3 = aJ(<p). 

As a check for this method when we choose the function J(0) to be constant (say, when J((f>) = 1), we can see 
that these equations represent the special case of a Minkowski universe, known to be a fixed point in the phase space 
describing the dynamics of the system. 

Another import check is to see if the case of an universe with an exponential expansion is recovered when we set 
W{<p) = const, as it is expected since 

H{4>)=W{4>) = H — > -=H — > a(t) = a(0)e Hot , (18) 

a 

but without imposing a form for a(t). We will try to find valid cosmological solutions using FOF for different values 
of n and different forms of W{4>) and J{<f>). 



IV. EXACT SOLUTIONS FOR THE INTERACTING DARK ENERGY AND DARK MATTER 

MODELS WITH FOF 

In this section, we search for cosmological viable solutions using the FOF. We show some examples of solutions 
giving different evolutions and couplings between the dark matter and dark energy. We intend to study how to obtain 
the cosmological solutions for different couplings and understand how well the method can describe the coupling 
scenario. As in the case of the phase space analysis we expect to find at least one cosmological viable solution 
presenting a period where the dark energy energy density becomes subdominant. In all cases I = 1. 

After finding the solutions by the method, the free parameters of each of them must be constrained by the ob- 
servational data in order to determine in which circunstances the solutions are viable. Some solutions required the 
introduction of parameters in order to have the right dimensions, that were also matched with the observations. 

In the first two examples we describe an accelerated universe, given W = const., and see how the FOF for different 
n and different J((f>) gives us different evolution and couplings between the dark energy and dark matter. This also 
serves as a test for the ansatz since we know that for that W(<f>), the scale factor must evolve exponentially. We are 
interested in verifying if the solutions obtained are cosmological viable by evaluating the cosmological parameters for 
each example. 



A. 1 st Example: 

Considering the case with n = 0,cr = 1 and W{(f>) = Hq, J{4>) = ~ q M , where C is a constant with dimension of 
energy (MeV), the Eqs. (|lHll4j) result in: 

F (4>) = -a + ° ~ R , V (4>) = -THq + — ± (0 - Mf , (19) 

p ap 2 2 

/ r \ 1/3 

4>(t) = M — Dc~ 3Hot , a(<)= (^) eH0t - (20) 

Such a universe develops a de Sitter expansion, therefore an accelerated expansion. We can also see this by 
calculating the acceleration parameter q = tia/a 2 . In this case q is constant and equals to 1. The interaction, 
described by F(<f>) and that also appears in V(^), accounts for the interactions that were introduced between the dark 
sectors and the one introduced by the FOF, as pointed out in [4lj . 

Here, the Hubble constant Hq and the a parameter play the role of a coupling constant between dark matter 
and dark energy through the relation (3p = — 9H ° . These quantities also redefine the mass of the dark matter as 
rriF = — 9H ° AI , making it tachyonic. 
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In order to understand the role of such an extra interaction, the cosmological model obtained and to confirm 
the accelerated expansion seen in the acceleration parameter, we evaluate the energy density, the pressure and the 
equation of state parameter for the scalar field, 

^DE(a) = 1 + £^p^6i PDE{a) = -^TH^ w DE (a) = o ° — ■ (21) 



We can also obtain these quantities for the dark matter, 



9ggg 2 _l 3 r rr2 
2 1 -"0 



n , . 2aM 1 . , 9H2 1 aM , , 9C# 2 1 

<Wa)=3^^, Plw (a) = -_--_, WDM ( a ) = _!____, (22 ) 



'0 

and for the interaction 



n f s 6 12 aM 1 Pint(a) , \ -. /OQ ^ 

n«(a) = ^ " s^f ^ = -7-7- • = - 1 ■ ( 23 ) 

We can readily see that i7 tot = fi^s + Qdm + ^int = 1 as given by the Friedmann equation 

Given the form of the solution, we need to determine aM and C by adjusting the solution to the observational 
data. We constrain some of the constants by fixing the density parameter of the dark matter today- a,Q = 1, the value 
of the scalar factor today- to be approximately 0.26, according to WMAP-7 data [47|. This fixes aM. 

We still need to constrain C. Since C 2 that appears in the energy density and in the equation of state cannot be 
negative so the scale factor is not imaginary, we can see from the expression for density parameter of dark energy that 
this will always be greater than one. Adjusting it to be as close to one as possible, so we do not have an over-density 
of dark energy, we were able to constrain C to a value of the order of 10~ 8 . 

In the left panel of Figure [T] we plot the density parameter of all the components. The interaction density is always 
negative and large. We can see that if we put this component together with any of the other components, how it is 
usual to do, this could influence a lot the final energy density of the components. 

We can see in the plot that the component that the "dark matter" dominates in the past, while dark energy 
dominates at late times. However, this does not mean that we observe deceleration in the past, as we already know 
it does not happen because of the acceleration parameter equal to 1. As we can see in the right panel of Figure Q] 
the component we are calling dark matter actually accelerates the universe since its equation of state asymptotes —1 
towards the present, instead of being equal to zero as expected for dust. This happens for two reasons: first, because 
the mass term of dark matter is cancelled by the first factor in F{<j>) in pTjl and only the kinetic part plays a role; 
second, because of the ansatz chosen, the scalar field will play a role in all the terms. So, the evolution of the "dark 
matter" component contains the interaction and behaves differently than the usual dark matter. 

For the dark energy component we can see that the equation of state presents an initial period where this fluid 
would not accelerate the universe, because it is larger than —1/3. So, this equation of state tends to a "matter" one 
in the past and then, around a = 0.2, it enters in the regime where it accelerates the universe going to —1 in the 
present and behaving like dark energy. This feature of the equation of state can be found in the equation of state of a 



Chaplygin gas 51| , where one component can mimic the behaviour of matter and dark energy and it is an alternative 
to the quintessence scenario of dark energy [f54| . 

However, as the dark energy component is not dominant during the period where it does not accelerate the universe, 
period in which the "dark matter" component that accelerates is dominant, the effective behaviour is of an accelerated 
universe during all times. This can be seen also be the acceleration parameter, in the right panel of Figure [TJ It is 
important to notice that this solution describes a flat universe (the sum of the density parameter of the components 
is equal to one) that is accelerating. 

Because this solution describes a universe that is only expanding in an accelerated way, this solution is not a good 
candidate for a viable cosmological solution for the late universe. The fact that we could only obtain an accelerated 
solution should not be used as a reason to discard it, since we imposed it by setting the Hubble parameter to be 
constant as an hypothesis. This solution could be a good one to describe a de Sitter inflationary period with two 
different components responsible for the acceleration in different times of the evolution. This solution is an attractor 
solution in the phase space analysis of dark energy (37| and appears naturally in our formalism, as would have been 
expected. 

We now show another solution, for a different n but with the same constant Hubble parameter. 
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FIG. 1: Density parameter f2 (left panel) and equation of state parameter w and acceleration parameter (right panel) for the 
model represented by eqs. I2H23I 



B. 2 nd Example: 



Considering the case with n = 1,<t = 1 and functions W((f>) = Ho and J{4>) = — sypfjrgjY^T] > wnere ^ is a 
parameter with dimension 4 in energy (MeV ), Equations (|1HI14I) result in 



m = (3Pt)V3 f 0(t) = ^aO^ 



1/3 

e Hot . (25) 



As in the previous example, this universe also develops a de Sitter accelerated expansion with acceleration parameter 
equal to 1. This can be shown by evaluating the acceleration parameter that is constant and equal to one, as can be 
seen in the right panel of Figure [U showing that the solution is always accelerated. 

The new coupling constant is /3f = Ci/3 and the mass of the dark matter is cancelled by the first term of the 
interaction F(<p). In this case dark matter is represented by massless fermions, i.e., they contribute only with a 
kinetic term. 

Another solution for <j)(t) is possible but this choice results in a contracting universe. As we are not interested in 
this type of solution, (although we find it interesting that the formalism also presents us this sort of situation) we 
discard it. Here, the influence of the dark matter is explicitly present in the scale factor. 

For this model the energy density, the pressure and the equation of state parameter are given by 

2 / p \ I l 3 —1 

^DE(a) = 1 + — [——) -, p DE (a) = --THq , w DE (a) = 5 , (26) 

31 \3H J [ln( 7 a)]3 2 i , 2 / p * 



3r \3H J [ln ( 7Q )]3 



where 7 3 = P 2 /daf3. 

We can also obtain these quantities for the dark matter and the interaction, 



n , . 2Ma 1 , . fPH2\ a 1 Ma 

it dm {a) = OTT2r ^ > PDM{a) - 



o tt2t^ 3 ' rum v."v I o I 4 q i 

ZHgTa 3 V 3 J [ln(7a)] 3 a 

udmM = -I - (— V — ^— r, (27) 

' \3HoJ Ma [ln(7a)] l 1 ' 



and for the interaction 



U int (a) = -—[—=-) --^T7TF^ = ' w mt (a) = -I. 28 

3r \3H J [ln( 7 a)] 3 3H a T a Pcrit 
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We can constrain the constants in this example, P, 7 and aM, like we did previously by fixing the density parameter 
of the dark matter today to be approximately 0.26, according to WMAP-7 data [47j, to fix aM. Setting a — 1 today, 
we can constrain the coefficient 7 to be of order of 10. As the density parameter of dark energy is always greater 
than one, the remaining constant, namely P, cannot be fixed using the density parameter of today, since P 2 cannot 
be negative so the scale factor is not negative. To constrain this constant we use the equation of state of dark energy 
that must be negative for all times, asymptotically approaching -1, as shown in the right panel of FigureO However, 
in order to u>de ot — 1 at do = 1, P 2 must be small. Hence, in order to impose this condition we fine-tune P to be of 
order of 10 -7 to give us a realistic cosmological evolution. 

In the left panel of Figure [2] we plot the density parameter of the component that we called dark matter, which 
dominates in the past, and of dark energy that dominates in the future. The density parameter of dark energy is 
always greater than one but it is compensated, as it is the dark matter one, by the interaction that it is always 
negative, giving a flat accelerated solution, with f2 tot = Qde + &dm + &int = 1- This interaction is strong and if 
included in any of the other components, it would greatly affect the evolution. 

As in the first example, this does not mean that we have an early period of deceleration. As we can see in the right 
panel of Figure [21 the equation of state of the "dark matter" is negative, tending to —1 through all the evolution. This 
component accelerates the universe, different than what is expected from dust. Again, because in this example we 
have a cancellation of the dark matter mass and because the evolution is mixed with the scalar field one, the effective 
behaviour is similar to the presence of a cosmological constant. 

Also as the first example, the dark energy component decelerates in the past, starting to accelerates around a = 0.15 
as uj reaches —1/3 and continuing accelerating until the present with the equation of state tending to — 1 . This is also 
similar to the Chaplygin gas, since it evolves like dust and dark energy for different times. 

The effective behavior is that we only observe acceleration, as seen in the acceleration parameter, since the "dark 
matter" component that accelerates the universe dominates during the period where the dark energy decelerates. 
Then dark energy dominates and continues to accelerate the universe. 

This solution, again, is always accelerated, as expected since the Hubble parameter is constant. It could describe 
only the period of an inflationary universe with two accelerating components ath different times, for example, and it 
is not a good description for the dark energy in the late universe. 



a 




FIG. 2: Density parameter Q. (left panel) and equation of state parameter w and acceleration parameter (right panel) for the 
model represented by Eqs. (|26H2"51l . 

We present now other models in the attempt to check realistic cosmological solutions for non-constant Hubble 
parameters. 



C. 3 rd Example: 



Considering the case with n = l,cr = 1 and functions W(4>) = + Ho and J(<f>) = — 1/C 2 , where C is a parameter 
with dimension of energy (MeV), Equations (| 1 HI 14[) result in 

iC 2 (j> f I + H 

F (4>) = \ '- + C x , (29) 

pa 
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Vy{4>) = {4a 2 (-fid + Mf + C 2 [4C 2 (3H 2 T - T 2 ) + 24H o 0(a/3Ci + C 2 T ~ aM) 
+120 2 (C 2 (3H 2 + T) + a {pd - M)) + 36H O C* 2 3 + 9C*V] } , 



B I tB 

4>(t) = -ff + -tanh — 



B = J9H 2 d 2 + 6(-dl3a - C 2 T + aM) , 



(30) 
(31) 



o(t) 



6C 4 , o / iB 
C0Sh 2C 



B 



1/3 



(32) 



In this case the mass of the fermions were not cancelled and the interaction has two effective coupling constants. 
The acceleration parameter is given by 



q{a) 



Sa 3 - 1 



(33) 



where S = —B/6C 4 . This acceleration parameter can present different behaviours depending on the value of 5. This 
factor will be determined by the observational constraints on the parameters. 
For this universe, the energy density for dark energy is: 



n DE (a) = I [C 2 {2Y 3H 2 ) + 2a(Pd M)] - 24^g^i 
i- BC- 



Sa 3 



Sa 3 - 1 



B 2 C 2 

3r 



Sa 3 - 1 



C 



Sa 3 - 1 



Sa 3 
Sa 3 



3B 2 T 

where we do not show the other parameters because they are too lengthy. 
We can also obtain these quantities for the dark matter, 



p- { [9tf 4 + 12(C 2 - l)H 2 T - 4r 2 C 2 ] + 4a 2 (M - pdf + !2aH 2 (M - (3d)} jj—^ 



(34) 



2aM 1 

C 2 T 6a 3 - 1 ' 
1 



UDM\a) 



6aM 



Vdm{o) 



B' 1 



1 



36C 4 5a 3 



, 4 (<k 3 -l) 



6aPd~9C 2 H 2 + B 2 d iSa 



QaPd - 9c z H£ + B Z G 
3 1 " 



Sa 3 



6a 3 



(35) 



and for the interaction, 



{Sa 3 ) ] - 



3C 2 r 5a 3 - 1 



6a/3d-9C"H£+B 2 C 4 



{Sa 3 



Sa 3 



Pint (a) 

Pcrit 



w int {a) 



-1 



(36) 



We have to fix the parameters /3Ci, aM and C using the observational data. From (|33p . depending on the value 
of (5 we can have periods of only acceleration and only deceleration or periods that combine both. We choose a S 
parameter that gives us an early period of deceleration followed by a period of acceleration, as it can be seen in the 
right panel of Figure [3] This is chosen because this is the expected behaviour for the late evolution of our universe 
(with the restriction that here it is only composed only by dark energy and dark matter). 

To help fixing the remaining parameters, we have adjusted the density parameter of the dark matter today to be 
approximately 0.26 and the density parameter of dark energy to be approximately 0.74, according to WMAP-7 data 
[47l | adjusted for only two fluids. This way we were able to fine-tune the combinations of the constants /3Ci, aM and 
C in (|3"4"|) and (JHSJ) so the cosmological parameters agree with the observations. We plotted the density parameter 
and equation of state. 

In the left panel of Figure [3] we can see that dark matter dominates in the past, reaching its today value as a = 1, 
followed by a late period of dark energy domination, that also has the right density value today. Like seen in the 
acceleration parameter in the right panel of Figure [31 where we have deceleration followed by acceleration, we can see 
from the equation of state that the dark matter component decelerates the universe during all periods while the dark 
energy one accelerates the universe during all times. The only exception is during the very early universe where the 
equation of state of dark matter, that is plotted normalized for comparison reasons, is negative. However, this phase 
is very short and our simplified model with only two components is not a good approximation for this early period. 
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The interaction density is again always negative. However, different from the other examples it is small compared to 
the others. In this case the interaction would not spoil the behaviour of other components. 

So, this model is a viable cosmological solution, since we have an early period of dark matter domination followed 
by a late period of dark energy acceleration. It describes a flat universe since fltot = ^de |!2dm + ^mt = 1) similar to 
the standard cold dark matter plus a quintessence field scenario. This was obtained by fine-tuning the free parameters 
of the solution. However, the early evolution is different from this standard scenario which may indicate that with 
this choice of parameters the interaction is important only in the early universe. This shows the power of the method 
to provide good candidate models to our universe. 




FIG. 3: Density parameter (left panel) and equation of state parameter w and acceleration parameter (right panel) for the 
model represented by Eqs. (|29l36fl . The equation of state of dark matter is normalized in order to make the graph easier to see 
in a way still holds its main features. 

As we can see in this example, the role of the parameters is decisive in the behaviour of the solution. The fact that 
the interaction gives large contributions in the first and second examples could be what causes the solution to not 
exhibit any other behaviour than acceleration, since the formalism gives a higher weight in the scalar field component, 
by making all the assumptions on H and a(t) as depending on cf>. This contribution is imposed by the FOF,as pointed 
in |4lj . making the interaction large when adjusted to the observational data and spoiling the cosmological solution. 
In this third example, we obtained a good cosmological solution and the contribution from the interaction was small. 



D. 4 th Example: 



We will further investigate good cosmological solutions given by the method. We present now an example that 
exhibits a viable cosmological evolution but with different features than the previous one. 

Considering the case with n = 3,er = 1 and functions W((f>) — and J(<f>) = — 4 a p(M~^F(0)) > wncrc -P is a 
parameter with dimension four in energy (MeV 4 ), the Eqs. (|lHITlj) result in 



^'^ (apt) 2 / 3 



e *r . (38) 



This is also a case presenting massless fermionic dark matter, and an interaction displaying a product of exponential 
and inverse powerdaw functions, with coupling constant /3p = C\[3. We can obtain from this solution the acceleration 
parameter 

r / 6 \ 1/3 

«^-Hw) (39) 
= 1 -2kky (40) 
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where £ = (2P 2 /a/3Ci) 1 / 3 . We can see from this expression that for 



t > 



(2T) 



3/2 



6P 



(41) 



the expansion is accelerated, given that q > 0. So, this solution presents a transition between decelerated and 
accelerated expansion as we can see in Figure [H what is expected by the observations. 

For this model the energy density, the pressure and the equation of state parameter for dark energy are given by: 



^DE(a) = 1 



1 



31n(£a 

For the dark matter and the interaction, 
8r 2 Ma ln(£a) 



Pde{o) 



and for the interaction 



Vdm{o) 



P 1 



1 



8r 2 ln(£a) 



P 1 



1 



8r 2 (ln(£a)) 5 



(1 - 31n(£a)) , ude{o) 



aM 

— u>dm [a) 
£a' 3 



(l-31n(£a)) 
(l + 31n(£a)) ' 



8r 2 aM ln(£a) ' 



(42) 



(43) 



8T 2 Malaga) Pmt (a) 



31n(£a) 3P 2 



Pcrit 



Wint(a) 



-1 



(44) 



The parameters we have to fix in this example are £, P and aM. First, we can fix the parameters by imposing 
that the transition from deceleration and acceleration happens in the near past, as expected for a good cosmological 
model, choosing £ ~ 2, 02. We also adjusted the density parameter of the dark matter today to be approximately 
0.26, fine-tuning the ratio aM/P. 

With that we can evaluate the density parameters, plotted in the left panel of Figure 2] The density parameter 
of dark energy is always greater than one and dominates through all evolution. The dark matter parameter is also 
plotted and is subdominant and the interaction is always negative and is small compared to the other components, 
except for very early times. 

Although dark energy dominates during all times, as it can be seen by the equation of state of dark energy in the 
right panel of Figure |4j the dark energy decelerates the expansion in early times. It changes from deceleration for 
acceleration when reaches uj = —1/3, in the same time as the acceleration parameter changes sign. This shows that 
the density and equation of state evolutions arc coherent with the evolution shown in the acceleration parameter. 

As already pointed out, this dark energy component that behaves like "matter" and dark energy is analogous to the 
Chaplygin gas and in this example it is responsible for the entire evolution of the universe mimicking its components. 

The so called dark matter component has a large negative value. This happens because of the interaction and 
because in this example, as in the first and in the second, the dark matter mass is cancelled. 

This solution represents a flat universe, as £Ide + !1dm + Qint = 1, presenting all the necessary requirements for a 
viable cosmological model and describes a universe where the same scalar field is responsible for the early inflationary 
acceleration and the late one. 

We notice that in this model the interaction given by the FOF is small, after adjusting the parameter to fit the 
observational data, and it represents a viable cosmologial solution. A possible reason for this is that the form one 
chooses for W((f>) and J{4>) leads to a coupling that can potentially spoil the cosmological solution, since only small 
couplings are cosmologically acceptable in models of interaction between dark matter and dark energy and since the 
FOF favours the energy density of the scalar field and large parameters enhance this component. 



V. CONCLUSIONS 



Interacting cosmological models have attracted much attention in the last few years, both as a theoretical laboratory 
as well as a describing phenomenology. However, due to the enormous difficulties in obtaining solutions for the 
Friedmann equations in such an approach, the theoretical effort has not gone too far. The best description of the 
possible solutions to these problems comes from the dynamical systems' analysis, which gives an overview of all the 
families of possible solutions. One could ask for specific exact solutions for a particular interaction, what is known to 
be very difficult, if not completely impossible, as is usual in the framework of General Relativity. 

The method described here has the purpose of simplifying the search for exact cosmological solutions for an inter- 
acting dark sector. Such a simplification takes place after a series of assumptions, that however do not weaken the 
value of the solutions since they can be suitable for some specific situations. 
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a 




FIG. 4: Density parameter $1 (left panel) and equation of state parameter w and acceleration parameter (right panel) for the 
model represented by eqs. 1371441 

We have obtained a series of solutions and discussed their properties. As in the phase space analysis method, we 
have obtained solutions where the dark energy were the only dominant component and two solutions where we describe 
both decelerated and accelerated periods. They are viable cosmological solutions since the deceleration period will 
allow for structure formation, and the accelerated one will account for the observed late accelerated expansion of the 
universe. 

However, although they present the same behaviour these solutions accomplish that in a different way. The third 
solution has a massive dark matter component with an equation of state that decelerates the universe, dominating in 
the past, and a dark energy component that accelerates in the late epoch. The fourth example has a non- massive dark 
matter component always subdominant and a dark energy component that plays the hole of a matter-like component , 
decelerating the universe in the past, and of dark energy, accelerating it in late times. This behaviour is similar to 
the Chaplygin gas, although it does not present a phantom epoch. 

The first and second solutions present an accelerating-only expansion. As a constant Hubble parameter has been 
imposed, they show that the method is robust. The interaction given by the method plus the imposed one are large 
in this examples. In the third and fourth examples, the ones that give viable cosmological solutions, the interaction 
is small, at least in most of the evolution. 

Hence, we have showed that the FOF is capable of giving viable cosmological solutions in the same way as the 
phase space analysis. It also gives the expected attractor solutions where only expansion is observed. 

The analysis of specific models, with well-motivated interaction functions and coupling constants, and the compar- 
ison with observational data will be the subject of future work. 
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Appendix: Direct Approach for Solving the System of Equations 

It is worthy to stress the difference between the FOF and the more direct approach to the solution of the cosmological 
equations. We shall do this by presenting the following example. 

We take F{4>) = (f> and choose the de Sitter solution (a/a = Ho). In this case the equation of motion for the scalar 
field can be written as 

£(4> + SHocj,) +V' = ^. (45) 
a 6 

Using Friedmann's equations, we can eliminate the potential V and solve for the scalar field <p, obtaining 

4>{t) = R\ + K 2 e- 3Hot + K 3 c-i Hat , (46) 
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where Ki, K2 and K3 are constants. Note that, since there are no restrictions on the sign of these constants, the 
requirement of invertibility for <f>(t) is not fulfilled for all the solutions. 

For a power-law scale factor a = Kt p , with K and p positive constants, we have for 4>{t) 



4>{t) = Y x + Y 2 



(In*) 5 



Y.lnt 



(47) 



where Y±, Y 2 and Y 3 are constants. This solution is clearly non-invertible. Thus the FOF and the direct method can 
be understood as complementary formalisms. 



[1] A. Riess et al., Astron. J. 116 (1998) 1009. 

[2] S. J. Perlmutter et al., Astrophys. J. 517 (1999) 565. 

[3] N. A. Bahcall, J. P. Ostriker, S. J. Perlmutter, and P. J. Steinhardt, Science, 284 (1999) 1481. 
[4] T. Padmanabhan, Phys. Rep. 380, (2003) 235. 

[5] M. Trodden, S.M. Carroll, TASI lectures: introduction to cosmology, astro-ph/0401547. 
[6] A. Linde, Inflation and string cosmology, hep-th/0503195. 
[7] J. Frieman, C. Hill, A. Stebbins, I. Waga, Phys. Rev. Lett. 75 (1995) 2077. 
[8] R. R. Caldwell, R. Dave, P. J. Steinhardt, Phys. Rev. Lett. 80 (1998) 1582. 
[9] P.J.E. Peebles, B. Ratra, Rev. Mod. Phys. 75 (2003) 559. 
[10] S.M. Carroll, Phys. Rev. Lett. 81 (1998) 3067. 

[11] I. Zlatev. b. Wang, P. J. Steinhardt, Phys. Rev. Lett. 82 (1999) 896. 
[12] K. Skenderis and P. K. Townsend, Phys. Lett. B468 (1999) 46. 

[13] O. DeWolfe, D.Z. Freedman, S. S. Gubser, A. Karch, Phys. Rev. D62 (2000) 046008. 
[14] K. Skenderis, P.K. Townsend, Phys. Rev. Lett. 96 (2006) 191301. 
[15] R. Rosenfeld, J. A. Frieman, JCAP 09 (2005) 003. 
[16] C. Wetterich, Astron. Astrophys. 301 (1995) 321. 

[17] A.V. Maccio, C. Quercellini, R. Mainini, L. Amendola, S.A. Bonometto, Phys. Rev. D 69 (2004) 123516; b. Amendola, M. 

Baldi, C. Wetterich, Phys. Rev. D 78 (2008) 023015. 
[18] M. Baldi, Mon. Not. Roy. Astron. Soc. 411 (2011) 1077; M. Baldi, Mon. Not. Roy. Astron. Soc. 420 (2012) 430. 
[19] K. Koyama, R. Maartens, Y.-S. Song, JCAP 0910 (2009) 017. 

[20] buis P. Chimento, Alejandro S. Jakubi, Diego Pavon, Winfried Zimdahl Phys. Rev. D67 (2003) 083513; Winfried Zimdahl, 
Diego Pavon Phys. Lett. B521 (2001) 133-138; W. Zimdahl, Int. J. Mod. Phys. D14 (2005) 2319. 

[21] B. Wang, Y. G. Gong, E. Abdalla, Phys. Lett. B624 (2005) 141; B. Wang, C. Y. bin, E. Abdalla, Phys. Lett. B637 
(2006) 357. 

[22] B. Wang, J. Zang, C. Y. bin, E. Abdalla, S. Micheletti, Nucl. Phys. B778 (2007) 69. 

[23] S. Das, P. S. Corasaniti, J. Khoury, Phys. Rev. D73 (2006) 083509. 

[24] G. Huey, B. Wandelt, Phys. Rev. D74 023519 (2006). 

[25] E. Abdalla, B. Wang, Phys. Lett. B651 (2007) 89. 

[26] C. Feng, B. Wang, Yungui Gong, Ru-Keng Su, JCAP 0709 (2007) 005. 

[27] Z. K. Guo, N. Ohta and S. Tsujikawa, Phys. Rev. D76 (2007) 023508. 

[28] E. Abdalla, R. Abramo, b. Sodre, B. Wang, Phys. Lett. B673 (2009) 107. 

[29] Z. Huang, B. Wang, E. Abdalla, R. Su JCAP 0605 (2006) 013; J. Shen, B. Wang, E. Abdalla, R. Su Phys. Lett. B609 
(2005) 200. 

[30] S. Micheletti, E. Abdalla, B. Wang Phys. Rev. D79 (2009) 123506. 
[31] S. Micheletti, JCAP 05 (2010) 009. 
[32] S. Micheletti, 1009.6198 [gr-qc]. 

[33] C. Feng, B. Wang, E. Abdalla, R. Su Phys. Lett. B665 (2008) 111. 
[34] B. Wang, C. bin, D. Pavon, E. Abdalla Phys. Lett. B662 (2008) 1. 
[35] D. Pavon, B. Wang Gen. Rel. Grav. 41 (2009) 1. 

[36] E. Abdalla, b. R. Abramo, J. C. C. de Souza Phys. Rev. D82 (2010) 023508. 

[37] E.Copeland, M. Sami, S. Tsujikawa, nt. J. Mod. Phys. D 15 (2006) 1753. 

[38] B. Gumjudpai, T. Naskar, M. Sami, S. Tsujikawa, JCAP 0506 (2005) 007. 

[39] b. Amendola, M. Quartin, S. Tsujikawa, I. Waga, Phys. Rev. D 74 (2006) 023525. 

[40] D. Bazeia, C.B. Gomes, b. bosano, R. Menezes, Phys. Lett. B633, (2006) 415. 

[41] D. Bazeia, b. bosano, J. J. Rodrigues, R. Rosenfeld, Eur. Phys. J. C55 (2008) 113. 

[42] R. Courant, "Differential and Integral Calculus", Vol. 1. John Wiley & Sons Inc. (1988). 

[43] A. Saa, E. Gunzig, b. Brenig, V. Faraoni, T.M. Rocha Filho, A. Figueiredo, Int. J. Theor. Phys. 40 (2001) 2295. 
[44] E. Elizalde, S. Nojiri, S. D. Odintsov, D. Saez-Gomez, V. Faraoni, Phys. Rev. D77 (2008) 106005. 
[45] Z. Zhu, M. Hu, J. S. Alcaniz, Y.-X. biu, Astron. & Astrophys. 483 (2008) 15-18. 

[46] S. bepe, J. borca, F. Pena and Y. Vasquez, Int. J. Mod. Phys. D 20, 2543 (2011) [arXiv: 1109.0559 [gr-qc]]. 



13 



[47] E. Komatsu et ai, Astrophys.J.Suppl. 192 (2011), 18. 
[48] Elisa G. M. Ferreira and E. Abdalla, in preparation. 

[49] A. Guth, "The Inflationary Universe: A Possible Solution To The Horizon And Flatness Problems," Phys. Rev. D 23, 347 
(1981). 

[50] R. H. Brandenberger, "Inflationary cosmology: Progress and problems," arXiv:hep-ph/9910410. 

[51] N. Bilic, G. B. Tupper and R. D. Viollier, Phys. Lett. B 535, 17 (2002) [astro-ph/0111325]; A. Y. .Kamenshchik, 
U. Moschella and V. Pasquier, Phys. Lett. B 511, 265 (2001) [gr-qc/0103004]; Phys. Rev. D 66, 043507 (2002) [gr- 
qc/0202064]. 

[52] Similar separation can be seen in [3] 

[53] Although we put a general function in the interacting term that will cover a great range of possible interactions, we are 
already restricting that the interaction in the fermionic sector depends on t ptp. If F((f>) — we recover the Yukawa coupling. 
[54] Chaplygin gas can also develop a period of phantom evolution as seen in [5l[ that we do not observe here. 



14 



